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DELAMINATION  GROWTH  ANALYSIS  IN  QUASI-ISOTROPIC  LAMINATES 
UNDER  LOADS  SIMULATING  LOW-VELOCITY  IMPACT 

K.  N.  Shivakumar*  and  W.  Elber** 

NASA  Langley  Research  Center 
Hampton,  Virginia  23665 

SUMMARY 

A  geometrically  nonlinear  finite-element  analysis  has  been  developed  to 
calculate  the  strain  energy  released  by  delaminating  plates  during  impact 
loading.  Only  the  first  mode  of  deformation,  which  is  equivalent  to  static 
deflection,  was  treated.  Both  the  impact  loading  and  delamination  in  the 
plate  were  assumed  to  be  axisymmetric.  The  strain  energy  release  rate  in 
peeling,  Gj,  and  shear  sliding,  Gjj,  modes  were  calculated  using  the  frac¬ 
ture  mechanics  crack  closure  technique.  Energy  release  rates  for  various 
delamination  sizes  and  locations  and  for  various  plate  configurations  and 
materials  were  compared.  The  analysis  indicated  that  shear  sliding  (Gjj)  was 
the  primary  mode  of  delaraination  growth.  The  analysis  also  indicated  that  the 
raidplane  (maximum  transverse  shear  stress  plane)  delamination  was  more 
critical  and  would  grow  first  before  any  other  delamination  of  the  same  size 
near  the  midplane  region.  The  delamination  growth  rate  was  higher  (neutrally 
stable)  for  a  low  toughness  (brittle)  matrix  and  slower  (stable)  for  high 
toughness  matrix.  The  energy  release  rate  in  the  peeling  mode,  Gj,  for  a 
near-surface  delamination  can  be  as  high  as  0.5GU,  and  can  contribute 
significantly  to  the  delamination  growth. 

INTRODUCTION 

Composite  laminates  made  with  high  strength  resin  matrices  are  brittle 
and  are  easily  damaged  by  transverse  impact  loads.  Delamination  damage 
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commonly  occurs  under  Impacts  by  hard  objects.  De lamina t tons  are  usually 
invisible  from  the  surface  but  still  cause  significant  reduction  in  compres¬ 
sive  residual  strength  of  the  laminate.  The  formation  and  growth  of  dclaroina- 
tions  in  laminated  plates  have  recently  been  investigated  under  static  loading 
(ref.  1).  That  study  concluded  that:  (1)  de lamina t tons  are  formed  in  a 
laminate  when  the  transverse  shear  around  the  contact  area  between  the 
impactor  and  the  plate  exceeds  a  critical  shear  strength;  (2)  these  delamina¬ 
tions  grow  if  the  shear  strain  energy  release  rate  exceeds  the  matrix  shear 
toughness.  This  report  contains  a  finite-element  formulation  to  predict  the 
growth  of  an  already  created  delaraination  in  a  circular  quasi— isotropic 
laminate  using  fracture  mechanics  criterLa.  The  model  delaminat tons  shown  in 
figure  1  were  used  to  assess  the  effects  of  delamination  size  and  location  on 
the  strain  energy  release  rate  and  delaminat ion  growth.  Although  backface 
damage  is  one  of  the  common  failure  modes  in  thin  plates  under  impacts*  that 
problem  is  not  addressed  here.  Instead,  possible  delaminations  in  thick 
plates  are  examined.  For  simplicity,  the  impact  event  is  assumed  to  be  equiv¬ 
alent  to  an  axisymmetric  plate  loaded  slowly  in  the  direction  of  impact.  This 
assumption  is  reasonable  for  quasi-isotropic  plates  impacted  at  the  center. 

Early  studies  at  NASA  Langley  (refs.  2  and  3)  of  low  velocity  impact  on 
8  to  32-ply  circular  laminates  have  indicated  that  laminates  may  undergo  large 
deflections,  which  makes  the  problem  geometrically  nonlinear.  The  objectives 
of  this  study  are:  (1)  to  develop  a  geometrically  nonlinear  axisymmetric 
finite-element  (F-E)  analysis  which  includes  the  crack  closure  technique  of 
calculating  strain  energy  release  rates,  (2)  to  identify  the  mode  (opening  or 
shear  sliding)  of  delamination  growth  due  to  low  velocity  impact,  (3)  to 
examine  the  influence  of  delaraination  size  and  location  and  plate  config¬ 
uration  on  strain  energy  release  rates,  and  (4)  to  examine  the  compression 
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Instability  of  the  near- surf ace  delaminated  layer  due  to  Impact.  Figure  l 
shows  delamlnat Ion  models  considered  In  this  study.  A  mldplane  delamination 
(fig.  1(a))  was  selected  to  Ldentlfy  the  mode  of  delamlnatlon  growth,  l.e., 
mode  1  or  IT  and/or  combined.  The  critical  location  of  the  delamlnatlon  was 
predicted  by  analyzing  a  plate  with  three  delamlnat  Ions ;  one  at  the  midplane 
and  the  other  two  located  at  equal  distances  on  each  side  of  the  mldplane 
(see  fig.  1(b)). 

Compressive  flexural  stresses  that  develop  during  the  Impact  of  a  thick, 
laminate  may  cause  local  instability  of  delaminated  surface  piles.  The  local 
buckling  of  the  surface  sublaminate  results  in  high  interlaminar  peeling 
stresses  in  addition  to  shear  stress  at  the  delamination  front.  The  mode  of 
delamlnatlon  growth  therefore  could  be  mode  I,  which  is  different  from  cases 
(a)  and  (b).  To  examine  the  energy  release  rates  related  to  local  insta¬ 
bility,  a  plate  with  a  delamlnatlon  nearer  to  the  top  surface  of  the  laminate, 
referred  to  as  a  near-surface  delamination,  was  analyzed.  Although  more  than 
one  delamination  may  be  present  or  created  during  the  impact,  for  simplicity 
they  are  neglected. 

A  geometrically  nonlinear  two-dimensional  finite-element  program  (ref.  4) 
was  extended  to  include  the  axlsymmetrlc  case.  The  analysis  includes  calcula¬ 
tion  of  strain  energy  release  rates  In  the  two  possible  fracture  modes,  namely 
the  peeling  (G^)  and  the  shear  sliding  modes  using  the  linear  elastic 

fracture  mechanics  crack  closure  technique  (ref.  5).  Strain  energy  release 
rates  for  various  delamlnatlon  sizes  and  locations  as  well  as  plate  configura¬ 
tions  were  calculated  and  compared. 
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of  the  delamlnat ton  1.9  represented  by  z^;  for  midpiane  delatntnat Ion  z^  ■  0. 
The  Impact  by  the  sphere  is  assumed  to  be  ac  the  plate  center  and  is  repre¬ 
sented  by  a  static  equivalent  force  P.  The  deformation  of  the  quasi¬ 
isotropic  circular  laminate  and  sublaminates  in  the  delaminated  region  are 
assumed  to  be  axisymiaetric .  The  radial  and  transverse  displacements  of  the 
plate  are  represented  by  u  and  w,  respectively.  The  contact  force  between 
the  sphere  and  the  plate  is  represented  by  an  elliptical  pressure  distribution 

c  ^ 

(ref.  6)  (see  fig.  3),  such  that  P  =  2 n  I  C  q(r)r  dr,  where  q(r)  is  the 

J0 

contact  pressure  at  r.  The  radius  Rc  of  the  pressure  distribution  depends 
on  the  impact  force  and  stiffness  and  curvature  properties  of  the  sphere  and 
the  plate  (see,  for  example,  ref.  7).  For  the  range  of  forces  and  the  sphere 
and  plate  considered  here,  the  value  of  Rc  varied  from  0  to  1  mm.  For 
simplicity,  a  nominal  radius  of  0.5  mm  was  assumed  throughout  the  analysis. 

ANALYSIS 

A  geometrically  nonlinear,  plane  stress  and  plane  strain  finite-element 
(P~E)  computer  code  GAMNAS  (ref.  A)  was  extended  to  include  axisymmetric 
analysis.  The  program  uses  a  four  noded  isoparametric  element.  The  method¬ 
ology  followed  in  the  analysis  is  given  in  reference  8.  A  finite-element 
model  for  a  midplane  delamination  is  shown  in  figure  3.  The  linear  elastic 
fracture  mechanics  (LEFM)  crack  closure  technique  (ref.  5)  was  employed  to 
calculate  scraln  energy  release  rates  in  mode  I  and  mode  II  at  the  delamina¬ 
tion  front.  A  separate  study  wa9  made  to  confirm  the  applicability  of  LEFM 
crack  closure  technique  for  calculating  strain  energy  release  rates  In  this 
geometrically  nonlinear  problem.  Figure  4(b)  shows  a  frec-body  diagram  of  the 
delaminat ion  front.  Displacements  in  the  r  and  z  directions  at  nodes  1 
and  2  behind  the  delamination  front  (node  0)  are,  respectively,  ,  U2,  w^ , 
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and  **2  (see  fig.  4(b)).  Delarainat ion  tip  nodal  forces  in  the  r  and  z 
directions  are  Fr  and  Fz,  respectively.  Then,  the  strain  energy  release 

fzIw2  “  wil 

rate  in  the  peeling  mode  (I)  is  GT  **  — 7— — -r —  and  in  the  shear  sliding 

I  4  nR  g  AR 

Fr I  u2  “  UJ 

mode  (II)  is  ^ — ,  in  which  R^  is  the  delamination  radius 
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and  AR  is  the  element  size  employed  near  the  delarainat ion  tip.  A  fine  F-E 
mesh  (<?R  “  0.05  mm  or  0.1  mm)  was  employed  around  the  delamination  front* 

Mesh  refinement  was  studied  to  select  a  value  of  AR  which  yielded  relative 
errors  less  than  1  percent.  Strain  energy  release  rates  for  various  delamina¬ 
tion  sizes,  locations,  and  plate  configurations  were  calculated  and  the 
results  are  presented  in  the  following  section. 

RESULTS  AND  DISCUSSION 

Strain  energy  release  rates  in  mode  I  and  mode  II  for  the  typical  delami¬ 
nation  configurations  shown  in  figure  1  were  calculated  using  the  linear 
elastic  fracture  mechanics  crack  closure  criteria.  The  plate  was  assumed  to 
be  loaded  statically  by  an  equivalent  impact  force  acting  over  a  radius  Rc 
(fig.  3).  Energy  release  rates  for  various  delarainat ion  sizes  and  location, 
and  plate  configurations  and  forces  were  calculated.  Elastic  material  proper¬ 
ties  used  in  the  analysis  are  given  In  table  I.  Results  are  discussed  in  the 
following  sections. 

A  comparison  of  calculated  Gj  and  GIf  for  the  delamination  configura¬ 
tions  shown  In  figure  1  indicated  that  for  delaminations  away  from  the  sur¬ 
faces  (figs.  1(a)  and  1(b))  Cj  is  very  small,  less  than  0.1  percent  of 
Cjj  based  on  a  def ormed-conf iguration  analysis.  Hence,  G^  results  are 
presented  only  for  near  surface  delansinations. 
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Midplanc  Delaminat ion 

Figure  5  shows  the  force  P  versus  the  plate  center  deflection  wc  for 
eight  delaminat ion  radii  (0  to  12.7  mm).  As  expected,  the  increase  in  delami— 
nation  radius  reduces  the  plate  flexural  stiffness,  hence  the  initial  slope  of 
the  force-def lection  curve  decreases. 

Figure  6  shows  the  variation  of  Cn  with  force  P  for  a  range  of 
de laminat ion  radii  from  0*635  to  6*35  ram*  An  examination  of  the  curves  for 
various  shows  that  for  a  constant  force  P,  decreases  with  an 

increase  in  delaraination  radius  R6  >  1.27  mm.  Therefore,  the  delaminat ion 
growth  stops  after  the  initial  growth.  It  is  noted  here  that  for  Rc  less 

D 

than  or  equal  to  Rc  +  h/2,  the  top  compressive  stresses  reduce  the  effect  of 
the  shear  stress  singularity  at  the  delaraination  front,  with  the  result  that 
calculated  Gn  values  were  small  (see  dashed  curve  «  0.635  mm). 

A  method  of  predicting  the  delaminat ten  growth  during  loading  is 
explained  in  figure  6.  If  a  plate  having  mode  II  fracture  toughness,  say, 

GIIc  «  1.0  kJ/ra2  and  R&  »  1.27  ( R&  >  Rc  +  h/2)  is  loaded  from  zero  by  a 

steel  sphere,  the  Gj^  increases  from  2ero  along  the  R^  -  1.27  curve 
(see  fig.  6).  At  P  *  0.87  kN,  Gjj  becomes  equal  to  CjIc  (1.0  kJ/ra2). 

Further  increase  in  the  force  P  causes  delaraination  to  grow  at  constant 
GII  “  GIIc*  The  locus  °f  Gn  is  shown  by  a  broken  line  in  figure  6.  The 
plot  of  R6  versus  critical  force  P  for  a  GUc  of  a  laminate  Is  hereafter 
referred  to  as  "delamination  growth  curve." 

Delaraination  growth  curves  for  three  materials  whose  elastic  stiffness 
properties  are  the  same  as  those  for  T300/520S  and  with  CUc  values  of  0.5, 

1.0,  and  1.5  kJ/m**  are  shown  In  figure  7*  Each  curve  gives  two  values  of 
^6  ^or  a  force  P.  The  smaller  value  corresponds  to  the  case  where 

the  delamination  lies  within  the  impact  pressure  region  and  the  larger  value 
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is  for  outside  the  region.  Tor  Us  x  the  compress l  ve  normal  stress  redures 

the  shear  stress  singularity  at  the  de  1  am  I  nnt  ( on  front,  bunco  larger  force  is 
required  for  ( he  delamination  to  heroine  C,,  er  it  leal.  Tills  force  could 
exceed  the  material  strength  In  transverse  shear,  tension,  and/or  crushing 
(see  re! .  1).  Therefore,  the  creation  and  Initial  growth  of  de laml nat Ions  In 

and  around  the  loaded  region  have  t o  he  governed  by  some  other  criteria  like 
Interlaminar  strength  criteria.  The  horizontal  dash-dot  line  shown  In  the 
figure  Is  the  boundary  (R  s  >  l.l  mm)  above  which  the  present  analysis  (LF.KM) 
is  valid.  Failure  ( n  the  region  below  the  broken  line  Is  controlled  by 
interlaminar  strength  criteria.  A  large  del  am  l  tint  ion  R^  >  1.27  nun  grows 
i90riotdntc.il  lv  with  the  force.  Furthermore,  the  rate  of  dpl.imtn.it  ion  growth 
(dR  v'dP)  Is  greater  for  low  values  of  Gjlc  than  for  high  values.  Low  values 
of  Ojj,  designate  low  toughness  material* 

The  effect  of  impact  force  on  energy  release  rate  Is  shown  In  figure  8. 
Results  for  21  A  mm  and  plate  sizes  12.7  *  2  (radius  x  thickness  In  milli¬ 
meters),  12.7  x  3,  and  2S.4  x  2  are  shown.  Calculated  Cjj  decreased  for 
bo tli  Increased  plate  thickness  (2  mm  to  3  ram  and  radius  12.7  mm)  and  plate 
radius  (12.7  mm  to  2S.4  mm)  from  the  reference  case  (12.7  x  2).  But,  Increas¬ 
ing  thickness  by  SO  percent  caused  a  uuch  greater  reduction  in  GtI  than  did 
doubling  the  plate's  radius. 

Figure  9  shows  Gn  versus  force  for  T300/5203,  C6000/HX205,  C600G/F1S5, 
and  7300/P l 700  composites.  Calculated  Gir*8  for  all  four  materials  arc 
practically  the  same,  even  though  their  stiffnesses  differ  (see  table  l). 
However,  a  de-lamination  in  a  brittle  composite  (like  T300/5200),  which  will 
have  lower  G*Ic,  will  become  critical  and  grow  at  lower  force  than  In  tougher 
composites. 
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Multiple  Delaminat Ions 

A  laminate  generally  develops  more  than  one  delamination  In  thickness 
direction  during  the  impact  (see,  for  example,  fig.  1(c)).  A  question  is: 
which  one  of  these  delaminations  is  critical  in  the  subsequent  impact  or 
transverse  loading?  To  predict  the  location  of  a  critical  delamination,  a 
plate  with  three  delaminat ions  was  analyzed.  One  of  the  delaminations  is  at 
the  midplane  (z§  -  0)  and  among  the  remaining  two,  one  is  assumed  to  be  above 
the  midplane  and  the  other  below  the  midplane.  Figure  10  shows  strain  energy 
release  rate  for  three  delaminat ions  (R^  -  4  mm)  located  at  -  0  and 
±0.5  mm.  The  results  indicate  that  Gu  for  the  midplane  delaraination  is 
higher  than  that  for  the  other  two.  Therefore,  a  delaaination  in  the  mid- 
plane,  which  is  a  plane  of  maximum  transverse  shear  stress,  is  most  critical 
and  grows  before  other  delaminations  in  that  region.  It  may  be  noted  that  the 
present  analysis  is  based  on  homogeneous  laminate  properties.  In  an  actual 
laminate,  maximum  shear  stress  location  depends  on  the  stacking  sequence  which 
is  a  function  of  layup  and  polar  angle  0  of  the  plate. 

Figure  11  shows  the  predicted  critical  forces  for  initiation  of  growth  of 
three  delaralnations.  The  material  fracture  toughness  GUc  is  assumed  to  be 
1.5  kJ/nT.  For  small  de laminat tons  (Rg  <  2.5  mm)  radii,  the  upper  delaaina¬ 
tion  (A)  grows  first  followed  by  B  and  C  delaminat ions.  However,  for 
R6  >  3  ma,  the  midp lane  delamination  grows  first  and  followed  by  A  and  C. 
These  curves  can  be  used  to  predict  delaminat ion  radii  for  a  given  force  by 
drawing  an  ordinate  perpendicular  to  the  abscissa  as  in  figure  11.  The  inter¬ 
section  of  the  ordinate  and  the  curves  gives  the  delamination  radii  at  the 
three  locations.  For  example,  an  impact  force  of  2.5  kK  causes  delamination 
of  radii  5.5,  4.4,  and  4  m  at  midplane,  z&  -  -0.5,  and  0.5  cm,  respectively. 
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Near-Surface  Dc  laiainat  ton 


Compressive  induced  instability  of  delaminated  laminas  located  near  the 
impacted  surface  in  a  thick  laminate  was  examined.  If  a  thin  debonded  sub- 
laminate  forms  near  the  top  surface  of  a  laminated  plate*  it  may  buckle  due  to 
the  compressive  bending  stresses  developed  by  the  transverse  force  on  the 
plate  (see  fig.  1(c)).  Large  peel  stresses  caused  by  local  buckling  would 
tend  to  produce  growth  by  mode  I  and  mode  II.  Figure  12  shows  energy  release 
rates  in  both  modes  for  delaminations  located  at  z&  -0.75  mm  at  a  force 
P  «  2.4  kN.  The  two  curves  Indicate  that,  due  to  local  buckling,  the  ratio 
of  Gl  to  Cir  couLd  be  as  high  as  0.5  (Rfi  «  4  mm).  Since  CIc  of  a 
material  Is  much  smaller  than  Gjjc,  the  delamt nation  growth  could  be  due  to 
peeling  alone*  Both  Gj  and  decrease  with  increase  in  R^  (R^  >  4  mm). 

Therefore,  the  delamination  growth  stops  after  an  initial  growth  for  a  given 
force.  For  small  values  of  R&  (R&  <  2  mm)  Gt  decreases  and  Gu 
increases.  Although  more  than  one  delamination  can  exist  in  a  laminated 
plate,  for  simplicity  and  to  narrow  attention  to  compression  buckling,  only  a 
near  surface  delaminat ion  was  considered  in  the  above  analysis. 

Fro©  figures  10  and  12  the  values  of  Gjj,  for  P  “  2.4  kN  and 

R5  •  4  mm,  and  for  delaminat ions  located  at  various  values  of  z g  are  given 
below. 


z  ^  mm 

-0.75 

-0.5 

0.0 

0.5 

Gn ,  kJ/m** 

0.21 

i 

l.l 

1.4 

1.0 

These  results  indicate  that  the  delaminat ions  near  the  esidplane  would  have 
reached  critical  Gn  and  grown  before  a  near-surface  delamination  becomes 
critical.  The  growth  of  near-surface  delamination  in  mode  II  is,  therefore. 


less  probable.  However,  the  existence  of  C£  mode,  even  though  small,  can 

cause  surface  delamlnatlon  growth  since  a  material's  CIc  is  far  smaller  than 
Is  GIIc. 

CONCLUSIONS 

An  axlsymraetric  finite— element  analysis  of  a  quasi-lsotroplc  circular 
laminate  under  a  static  equivalent  Impact  force  led  to  the  following 
conclusions; 

1.  The  primary  mode  of  delamlnatlon  growth  In  the  interior,  away  from 
top  and  bottom  surfaces,  of  a  circular,  quasi-lsotroplc  laminate  subjected  to 
transverse  Impact  Is  by  Gn  (Gm  was  not  addressed). 

2.  For  a  constant  force,  GtI  decreases  with  Increasing  radius  of  the 
delamlnatlon  when  the  delamlnatlon  radius  is  greater  than  the  loaded  region. 

3.  Small  delaminations  (smaller  than  or  nearly  equal  to  the  Impact  pres¬ 
sure  area)  have  very  low  C^'s  due  to  a  reduced  shear  stress  singularity. 
Therefore,  Initial  Interlaminar  failures  may  be  governed  by  the  interlaminar 
strength  criteria* 

4.  Delamlnatlon  growth  rate  Is  rapid  (neutrally  stable)  for  low  fracture 
toughness  (brittle)  mater  -'.s,  and  It  Is  slow  and  stable  for  high  toughness 
materials. 

3.  A  midp lane  delamlnatlon,  with  a  radius  larger  than  the  impact  pres¬ 
sure  radius.  Is  more  critical  and  grows  faster  than  other  delarai nat ions  of  the 
same  size  but  located  in  any  other  plane. 

6.  If  a  thin  dobonded  sublaminate  forms  near  the  Impact  surface  of  a 
thick  laminate,  the  sublaminate  may  buckle  due  to  bending  compressive 
stresses.  The  peel  mode  energy  release  rate  Gt  can  be  as  high  as  O.SGjj. 
Therefore,  a  near-surface  delamlnatlon  in  a  material  with  low  GIc  may  grow 
in  the  peeling  mode. 
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Table  1-  Elastic  properties  of  materials. 


Mater 


T300/5208 

C6000/HX205 

54.00 

48.33 

14.30 

9.17 

3.50 

4.10 

0.28 

0.27 

0.31 

0.31 

Lai 


C6000/F185 


T300/P1700 


44.35 

6.79 

2.70 

0.27 

0.32 


37.56 

6.64 

2.62 

0.38 

0.32 


DELAMINATION  RADIUS,  Re,  mm 


force  versus 


dclanination,  T300/5 


T300/5208 


ORIGINAL  PAC!  T3 
OF  POOR  QUALITY 
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Effect  of  force  on  for  various  materials  for  a  mid-plane  dclamination  plate, 

(a  -  12.7  rism,  h  ®  1.0  mm,  R  »  3.8  nm) 


in  energy  release  rates  for  three  delnsni nation  locations. 
»  12.7  tnn,  h  *  2.0  ran,  Rf  ^.0  bw,  T300/5208) 


mm 
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Fig.  11  Effc?ct  of  force  on  radius  of  dcinmination  for  multiple  delaminations, 

(a  -  12.7  ra,  h  *  2.0  am,  CTf  »  1.5  kJ/ta2,  T300/5208) 

lie 


energy  release  rates  for  near-surface  delaminations  for 
2.7  mm,  h  *  2  mm,  z  *  -  0.75  mm,  T300/5208) 


